Limits to error correction in quantum chaos 
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We study the correction of errors that have accumulated in an entangled state of spins as a result 
of unknown local variations in the Zeeman energy (B) and spin-spin interaction energy (J). A 
non-degenerate code with error rate k can recover the original state with high fidelity within a time 
ta ~ fiK 1 ' 2 /max(5, J) — independent of the number of encoded qubits. Whether the Hamiltonian is 
chaotic or not does not affect this time scale, but it does affect the complexity of the error-correcting 
code. 
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In classical mechanics, chaos severely limits the op- 
eration of a reversible computer Any uncertainty 
in the initial conditions is magnified exponentially by 
chaotic dynamics, rendering the outcome of the com- 
putation unpredictable. This is why practical compu- 
tational schemes are irreversible: Dissipation suppresses 
chaos and makes the computation robust to errors 0. 
A quantum computer does not have this option; It relies 
on the reversible unitary evolution of entangled quantum 
mechanical states, which does not tolerate dissipation ||. 
This invites the question [jU^] , what limitations quantum 
chaos might pose on quantum computing. 

To answer this question one needs to consider the pos- 
sibilities and restrictions of quantum error correction || . 
Errors can occur due to interaction with the environment 
(errors of decoherence) and due to uncertainty in the uni- 
tary evolution (unitary errors). The original state can be 
recovered reliably if the errors involve at most a fraction 
k < 0.1 of the total number of qubits. The correspond- 
ing maximal time during which errors may be allowed 
to accumulate (the recovery time £r) is easy to find if 
different qubits are affected independently. That may be 
a reasonable assumption for certain mechanisms of de- 
coherence and also for unitary errors resulting from an 
uncertain single-particle Hamiltonian. Uncertainties in 
the interactions among the qubits pose a more complex 
problem [Q. 

Georgeot and Shepelyansky Q] studied this problem for 
a model Hamiltonian of N interacting spins that exhibits 
a transition from regular dynamics (nearly isolated spins) 
to chaotic dynamics (strongly coupled spins). They con- 
cluded for the chaotic regime that % goes to zero oc 1/N 
for large N, but their analysis did not incorporate the 
optimal error-correcting procedure. We assume a good 
(non-degenerate) error-correcting code and obtain a re- 
covery time of the order of the inverse energy uncertainty 
per spin — irrespective of the number of encoded qubits. 
By considering both phase-shift and spin-flip errors we 
find that £r is insensitive to whether the Hamiltonian 
is chaotic or not. (The authors of Refs. |||f| arrived at 
the opposite conclusion that <r increases strongly when 
chaos is suppressed, but they took only spin-flip errors 



into account.) The absence of chaos can be used to re- 
duce the complexity of the code, in that a classical error- 
correcting code suffices for the majority of the errors in 
the regime of regular dynamics. 

The Hamiltonian H under consideration describes N 
coupled spins er„ = (ct^, a^,a^) on a lattice in a magnetic 
field B n = {B%,Bl,Bl), 



(1) 



A spin n interacts with d neigboring spins m via the 
matrix 3 nm . The spin could be a nuclear spin or the 
spin of an electron confined to a quantum dot, in the 
context of solid-state based proposals for quantum com- 
puting ||s|-|l0||. More generally, the spin could be a rep- 
resentation of a two-level system (for example, in the 
context of the ion-trap quantum computer ]TT[]). We 
assume that B„ and 3 nrn fluctuate independently from 
site to site, with zero mean and variance |B„| 2 = B 2 

and J2af3 (Jnm) 2 = J 2 (provided m is one of the d 
neighbors of n, otherwise J nm = 0). We denote by 
U = (B 2 + 2dJ 2 ) 1 / 2 the root-mean-square energy un- 
certainty per spin. 

A state ipo evolves in time according to ifj(t) = e~ %Ht i\)§ 
(setting Ti = 1). We assume that we do not know the 
parameters of the Hamiltonian, and use quantum error 
correction to recover ipo from i/j(t) [Q. Let ipo lie in the 
code space of a non-degenerate error-correcting code [0 . 
The code space is a 2 M dimensional subspace of the full 
2 N dimensional Hilbert space, such that 



•< fc W) = 0, l<k<2K, 



(2) 



for any two (possibly identical) states ipo, i/j' in the code 
space and any product of up to 2K Pauli matrices a" 
(acting on different spins n\, ri2, ■ ■ ■)■ The number M is 
the number of qubits encoded in N spins. The number K 
is the number of errors that the code can correct, where 
the application of a^,a^, or to any of the N spins 
counts as one error. The ratio M/N = p is the bit rate 
of the code and the ratio K/N ~ k the error rate. 
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Error correction is successful if tp(t) lies in the error 
space of tpo, which is the subspace spanned by the state 
ipo and the states derived from ip by making up to K 
errors. The operator V projects onto the error space. 



Explicitly: V = P P , with 



= E 

{n,a} 



(3) 



The symbol £{n a } indicates a summation over the n^s 
and aj's, with the restriction that the indices ni,ri2... 
should be all distinct. (The indices a±, a.2 ■ • ■ need not be 
distinct.) The norm 



F (t) = \vm\ 2 = (H 



AHt 



-iHt 



|V>o> 



(4) 



of the projected state is the probability of successful er- 
ror correction after a time t. It is the "fidelity" of the 
recovered state |(| . The recovery time £r can be defined 
as the time at which the fidelity has dropped from 1 to 
1/2. 

We assume that the error-correcting code is "good", 
meaning that p and n tend to a non-zero value as N — > oo. 
Good quantum-error correcting codes exist, but their 
construction for large N is a complex problem [^3| [L8j. 
Our strategy will be to derive a lower bound to F and 
£r that does not use any properties of the code beyond 
the non-degeneracy condition (^), so that we can avoid 
an explicit construction. An alternative approach would 
be to abandon the requirement of a "good" code, and 
keep the number M of encoded qubits fixed as the to- 
tal number of spins N goes to infinity. One can then 
use the technique of concatenation to construct codes 
that are safe for a large number of errors at the expense 
of a vanishingly small bit rate p. (See Ref. pl| for such 
a calculation in the case M — 1.) 

Our first step is to decompose the evolution operator 
e iHt _ Xk into operators Xk that create k errors. 

For k <C N and t -C 1/U we may approximate 



Xk — Xi 



(it) 



k-q/2 



asi 



k 



A ^711712 



TQ,_lQ, p>a q +l 
J -^n q+1 



B ak 



Xq = exp 



■l* 2 E 



n^m a.,/3 



(5) 
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The approximation consists in neglecting terms in the 
exponent of order k(Ut) 2 and iV(£/£) 4 , relative to the 
terms retained of order N(Ut) 2 . We may write Xq w 
cxp[— i7V([/i) 2 ], neglecting fluctuations in the exponent 

that are smaller by a factor 1/y/N. 

We next substitute the decomposition of e lHt in error 



operators into the fidelity (||), 



K N 

F (*) = E E ww>, 

p=0 k,k'=0 



(7) 



where we have abbreviated (• ■ ■) = (?/>o| ■ 1 1 IV'o)- To sim- 
plify this expression, we take the average over the random 
variations in the B n 's and Jnm's. (We will show later that 
statistical fluctuations around the average are insignifi- 
cant.) Only the terms with k = k' contribute to the 
average. The terms with p + k < 2K can be simplified 
further, since they contain at most 2K Pauli matrices. 
In view of Eq. (g), these expectation values vanish un- 
less the product of Pauli matrices reduces to a c-number, 
which requires p — k. Hence the average fidelity can be 
written as F — F\ + F-i , with 



Fi=J2 E ^\(X P a%l...a%W, (8) 

p=0 {n,a} P ' 
K N 

^=E E E ^\(Xk<r%...<r%)\ 2 . (9) 



p=0 k=2K +l-p {n,a} 
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FIG. 1. Solid curves: time dependence of the lower bound 
Fi to the ensemble-averaged fidelity, calculated from Eq. 
for error rate k = 0.01 and three values of N. (We took 
B 2 = 2d J 2 , so that the root-mean-squared energy uncertainty 
per spin U is equally divided between Zeeman energy and in- 
teraction energy.) The dashed curve shows (for N = 10 4 ) the 
squared overlap Xq = exp[— N(Ut) 2 ] between initial and final 
state. In all these curves the number M of encoded qubits is 
a fixed fraction p of the total number of spins N. 

The expectation values in F\ are evaluated by substi- 
tuting Eq. (^]) and extracting the terms that reduce to a 
c-number, 



=0 9=0,2,4 



(p-q)l (q/2)\ 



(10) 



For K > 1 we may approximate e~ x x k jk\ « 
exp[— (k — x) 2 /2x] and replace the sums in Eq. 
( |To| ) by integrals. The result is 



F x = \ + ±erf [(t R - t)/At 

*H = 



(Ha) 



B 2 + AdJ 2 



jl^B 2 + 8dJ 2 , , . 

' At = Vm b 2 + mj 2 ■ (llb) 
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with erf (a;) = 27r~ 1 / 2 f Q x e~ y2 dy the error function. Cor- 
rections to £r and At arising from the approximations 
made in Eqs. (pj) and (^) are smaller by a factor k. The 
expectation values in F% depend specifically on ipo, hence 
on the way in which M qubits are encoded in N spins. 
Since F2 > we have a lower bound F > F\ on the 
fidelity that is code independent within the class of non- 
degenerate error-correcting codes. 

For N — > 00 the time dependence of Eq. ( jO ) ap- 
proaches the step function 8(tR — t). The threshold t-R is 
independent of N, while the width At of the transition 
vanishes as iV -1 / 2 (solid curves in Fig. Q). These are 
results for the ensemble-averaged fidelity, but since the 
variance is bounded by < vaxF < F(l — F) the fluc- 
tuations are insignificant except in the narrow transition 
region. The step-function behavior of the fidelity also 
implies that the positive code-dependent term F 2 that 
we have not included in Fig. [l] satisfies limjv— »oo F2 —* 
for t < t R (since F x + F2 < 1 and F x -> 1 for t < i R ). 
Any code-dependence of the fidelity can therefore only 
appear for times greater than t-R. 

The independence of the recovery time on the number 
M of encoded qubits disagrees with Refs. 0,b[. These 
authors calculated the squared overlap | {"00 1 ^ (*) ) 1 2 ~ -^0 
of the time-dependent state with the original state, and 
argued that the original state would be effectively lost 
once this overlap is <C 1. However, the original state can 
be recovered even when this overlap has become expo- 
nentially small, if a good error-correcting code is used 
(compare dashed and solid curves in Fig. ||). The recov- 
ery time is increased by a factor V kM, with an overhead 
of 1/p spins per encoded qubit. 

We find that tR at a given U is insensitive to the rel- 
ative magnitude of B and J, and hence insensitive to 
whether the Hamiltonian is chaotic or not. This con- 
clusion may seem surprising in view of the fact that the 
eigenstates are completely different in the chaotic and 
regular regimes Q: For J < B/N the eigenstates of the 
total Hamiltonian H are a superposition of a small num- 
ber of eigenstates of the non-interacting part ^2 n B„ • cr„. 
This number (known as the participation ratio) increases 
with increasing J, and when J « B it becomes of the 
same order as the dimension 2 N of the entire Hilbert 
space. (See Ref. |20) for a description of the onset of 
quantum chaos in systems with random two-body inter- 
actions.) As we will now discuss, the reason that a small 
participation ratio does not improve the fidelity is that 
it counts spin- flip errors but not phase-shift errors. For 
the same reason, suppression of chaos does help to reduce 
the complexity of the error-correcting code. 

The three Pauli matrices correspond to three types of 
errors: spin flips (a x ), phase shifts (a z ), and a combina- 
tion of the two (a v = ia x a z ). The complexity of the code 
is reduced substantially if there is only one type of error 
to correct. (One can then use a code for classical bits, 
such as the Hamming code Q.) Suppose that we seek to 
suppress spin-flip errors, of either type a x or cr v . To this 



end we impose on the spins a known uniform magnetic 
field in the z-direction, with Zeeman energy Bq that is 
large compared to the magnitude U of the random en- 
ergy variations. The new Hamiltonian is H + Hq, with H 
given by Eq. (|j) and H = B J2 n °~n- Since B is known 
we can undo the evolution of a state due to Hq by apply- 
ing the operator e lHat — Y\ n (cos B$t + ia^ s'mBot). Any 
remaining deviation of i/j(t) from i/;q has to be dealt with 
by the error-correcting code, with projection operator V . 
The fidelity of the corrected state is F(t) = \PG(t)ip \ 2 , 
where the evolution operator G is defined by 



G(t) 



e iH t e ~i(H+H )t 



Texp 



1 I H(t')dt' 
la 



(12) 



(The notation T indicates time ordering of the operators 
H(t) = e iHot He- iHot .) 

For Bot > 1 we may replace H (t) by its time average 
over the interval (t,t + 1/Bq). The terms containing a 
single o~ x or o~ y average out to zero and we are left with 

G (t) = e"*W-), =J2 B n< + E J nm«, 

n n^m 

H±=E 5( J nm + J V Jn){« + «)■ (13) 



J y f n , so that the mixed 



(We have assumed 
terms <Jn°~m cancel.) The time dependence of the fi- 
delity is again given by Eq. (|j"T|), with B 2 = (B* L ) 2 and 

J 2 = (Jn z m) 2 + Wn?k + Jnl) 2 - The recovery time t R de- 
pends only weakly on the ratio J/B. The relative num- 
ber of phase-shift and spin-flip errors, however, depends 
strongly on this ratio. Indeed, if one would use a code 
that corrects up to Ku errors from a z and up to K± errors 
from a x or <r v , then the maximal £r (at fixed Kii + K±) 
is reached for K±/K^ = 4dJ 2 /B 2 . For J <C B one has 
K± -C K\\, so that a classical error-correcting code suf- 
fices for the majority of errors. 

Before concluding we briefly consider the case that 
the parameters B„ and 3 nm in the Hamiltonian are 
not only unknown but also time dependent. The re- 
sult (|10| ) still holds if we replace (Bt) 2 by the correlator 
b(t) = L dt' f Q dt" B„(i') • B„(t"), and similarly replace 

(Jt) 2 by j(t) = !ldt'f Q dt" E a0 J™(tV™(t")- For a 
short-time correlation one has bit) = bo\t\, j(t) = jo\t\. 
This leads for K > 1 to the fidelity 



±erf[(i r 



At 



t)/At], 



\/bo + 8dj 
N (b + 4^ ) 3 / 2 • 



(14a) 
(14b) 



The recovery time now depends linearly on the error rate 
k, but it remains iV-independent. The next step to- 
wards fault-tolerant computing would be to include in 
the Hamiltonian a part with a known time dependence, 
representing the logical gates. We leave that for a future 
investigation. 
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In conclusion, we have derived a code-independent 
lower bound for the fidelity F of a state that has been 
recovered after a unitary evolution for a time t in an un- 
known random magnetic field and spin-spin interaction. 
For a large system the transition from F = 1 to F = oc- 
curs abruptly at a time £r that is independent of the total 
number of spins N and the number of encoded qubits M. 
The magnitude of ts. is set by the inverse energy uncer- 
tainty per spin, regardless of whether the spins are nearly 
isolated or strongly coupled. The suppression of chaos 
that occurs when the spins are decoupled does not im- 
prove the fidelity, because of the persistence of phase-shift 
errors. Spin-flip errors can be suppressed, and this helps 
to reduce the complexity of the error-correcting code. 

In this work we have concentrated on the recovery from 
unitary errors. One might question whether suppression 
of quantum chaos improves the fidelity for recovery from 
errors of decoherence, in particular in view of the "hyper- 
sensitivity to perturbation" observed in computer simu- 
lations of systems with a chaotic dynamics pl| , ^2| . This 
question presents itself as an interesting topic for future 
research. 
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